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Abstract 

It is pointed out that an exactly solvable permutation operator, viewed as the quantization of cyclic shifts, 
is useful in constructing a basis in which to study the quantum baker's map, a paradigm system of quantum 
chaos. In the basis of this operator the eigenfunctions of the quantum baker's map are compressed by factors 
of around five or more. We show explicitly its connection to an operator that is closely related to the usual 
quantum baker's map. This permutation operator has interesting connections to the art of shuffling cards 
as well as to the quantum factoring algorithm of Shor via the quantum order finding one. Hence we point 
out that this well-known quantum algorithm makes crucial use of a quantum chaotic operator, or at least 
one that is close to the quantization of the left-shift, a closeness that we also explore quantitatively. 
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A textbook example of a simple fully chaotic system is provided by the model of a baker mixing 
dough, the baker's map. The classical baker's map 1], T, is the area preserving transformation 
of the unit square [0, 1) x [0, 1) onto itself, which takes a phase space point (q,p) to (q',p') where 
(q 1 = 2q, p' = p/2) if < q < 1/2 and (q' = 2q - 1, p' = (p + l)/2) if 1/2 < g < 1. The stretching 
along the horizontal q direction by a factor of two is compensated exactly by a compression in the 
vertical p direction. The repeated action of T on the unit square leaves the phase space mixed, 
this is well known to be a fully chaotic system that in a mathematically precise sense is as random 

Q 

as a coin toss |2j. The area-preserving property makes this map a model of chaotic two-degree of 
freedom Hamiltonian systems, and the Lyapunov exponent is log(2) per iteration. 

As the classical baker's map is exactly solvable in many ways, including an explicit prescription 
for finding periodic orbits of any period, its quantization was sought as a simple model of quantum 
chaos. The baker's map as many nice features, including 

simplicity, that make it ideal for this purpose and has been used extensively in studies of quantum 
chaos and semiclassical methods. It has also been experimentally implemented recently using NMR 
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. The quantum baker's map, in the position representation, that we use here is: 

\ G N/2 



(G N )W> = -L exp[-2vri(m + a)(n + (3)/N\. (2) 



where 



We require that /V be an even integer: Saraceno [f| imposed anti-periodic boundary conditions 
(a = (3 = 1/2) that we use. In this case we drop the superscripts indicating these phases. The 
Hilbert space is finite dimensional, the dimensionality N being the scaled inverse Planck constant 
[N = 1/h), where we have used that the phase-space area is unity. The position and momentum 
states are denoted as \qn) and |p m ), where m, n — 0, • • • , N — 1 and the transformation function 
between these bases is the finite Fourier transform G/v given above. 

The choice of anti-periodic boundary conditions fully preserves parity symmetry, here called R, 
which is such that R\q n ) = \qN-n-i)- Time-reversal symmetry is also present and implies in the 
context of the quantum baker's map that an overall phase can be chosen such that the momentum 
and position representations are complex conjugates: Gn4> = <jf, if </> is an eigenstate in the 
position basis. B is an unitary matrix, whose repeated application is the quantum version of the 
full left-shift of classical chaos. There is a semiclassical trace formula, which, based on the unstable 
periodic orbits, approximates eigenvalues 



Despite the simplicity of the quantum baker's map, its solution in terms of exact spectra con- 
tinues to be elusive. Recently we showed Q that for N that are powers of two, it is possible 
to write approximate analytic formulae for certain class of states. In particular the Thue-Morse 
sequence ({1,— 1,— 1, 1,-1, 1, 1,-1,...}, the n-th term is the parity of n when expressed in the 
binary, counting n from zero) Q] and its Fourier transform P] determined to a large extent a 
class of states we called "Thue-Morse states". Similar expressions were also found for families of 
strongly scarred states. Despite having simple, if approximate, analytic formulae these states were 
found to be multifractals. Thus we went some way in solving a quantum chaotic system that is 
nearly generic. A crucial tool used was the Walsh-Hadamard transform 10]. That is, if <j) is an 
eigenstate we studied Hk4>, where Hk = <8> H, a K-iold tensor product of the Hadamard matrix 
H = ((1, 1), (1, -l))/\/2, where 2 K = N. 

We wish to now address the case of general N and arrive at a counterpart of the Walsh- 
Hadamard transform that will simplify the states of the quantum baker's map. We show that a 
simple operator, the shift operator, that is exactly solvable, acts as a good zeroth order operator 
for the quantum baker's map. Therefore its eigenstates form a basis in which the eigenstates of the 
quantum baker's map appear simple. We study this operator's action in phase space, and show 
how to build a quantum baker's map around this operator. This "new" quantum baker will then 
turn out to be very close to the "usual" quantum baker's map in Eq. (Q). 

The shift operator S, by definition, acts on the position basis as S\q n ) = \q2n) or \q2n-N+i) 
depending on if n < N/2 or otherwise. We notice that S is "almost" B, only there is no momentum 
cut-off, as (p m \B\q n ) = V2(p m \Q2n) for n and m both < N/2 — 1. In fact S is a generalization of 
what was proposed as the quantum baker's map by Penrose for the case when N = 2 K . In this 
case if the position state \q n ) is denoted in terms of the binary expansion of n = o,k~\(ik-2 • • ■ 
then S\aK-idK-2 • • • clq) = \clk—2Q>K-3 • • • ao&KT-l)- It is easy to see that S commutes with the 
parity operator R. However S does not respect the usual time-reversal symmetry, relevant to the 
baker's map, namely GJ^S*Gn / S^ 1 . It does respect a "restricted" time-reversal symmetry 
in the case when iV = 2 , as b^ 1 S* b = S^ 1 , where b is the bit reversal operator defined as 
b\ a K-i a K~2 ' ' ' a o) = |°o a i " " " a K-2 a K-i) ■ It is useful to rewrite the action of S on the position 
basis (written simply as \n)) as 

S\n) = |2nmod(iV- 1)), (3) 

with the caveat that S\N — 1) = \N — 1), rather than |0). This is not crucial as it affects only an 
one-dimensional invariant subspace. 
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We point to two apparently unrelated contexts in which 5 has already appeared. Firstly S 
is closely related to the perfect "riffle-shuffle" (l^ used to randomize a deck of cards, to be more 
precise the "out-shuffle". If for instance N = 8 cards were in a deck, it is split into two exact halves 
and the cards are then interleaved. If the cards were numbered 0,1,2,3,4,5,6,7, the out-shuffle 
brings it to 0, 4, 1, 5, 2, 6, 3, 7, which is easily verified to be the action S~ l . The deterministic chaos 
of this shuffling process forms the basis of certain card tricks. The perfect shuffle returns the deck 
to its original state after a few shuffles, we will see below that this is the "quantum period function" 
relevant to S. 

Secondly, a generalization of S, where the factor 2 is replaced by any integer (coprime to 
N—l) is precisely the operator whose "phase estimation" leads to the solution of the order-finding 



problem 



13j |. The multiplicative order of 2 modulo N — 1 is the smallest integer r such that 



2 r = lmod(A r — 1), which is the quantum period again as S r = 1. We are guaranteed that such 
a number exists as Euler's generalization of Fermat's little theorem implies that (f>(N — 1) is such 
that 2*( Ar - 1 ) = 1 mod(iV - 1), thus r is either (f)(N - 1) or is a divisor of it (here 4>(n) is the Euler 
totient function, being the number of positive integers less than n and coprime to it). Finding the 
multiplicative order is the route of the quantum factoring algorithm of Shor. Thus it is interesting 
that this well-known quantum algorithm makes critical use of an operator that could be thought 
of as a quantization of the fully chaotic left-shift, or at least nearly, as explained below. 

That the classical limit of the unitary operator S is not the baker's map is made clear by studying 
its action on coherent states. The structure of S in the position basis is that of a permutation, and 
its action on the momentum basis is found easily: 

l m '\ ?\m\ - 1 ~ sm br(m' + l/2)/jVj + (-l) m+1 cos [Am' + l/2)/iV] 
{ m\b\m)- N sin[7r(m-2m'-l/2)/iV] {) 

Thus the momentum representation is also real. More importantly for a given initial momentum m, 
there are two momentum values around which the final state is spread, namely [m/2] or [m/2]±iV/2. 
Thus the action of S on coherent states would be roughly a combination of its actions on position 
and momentum states, and therefore splits an initial state while performing appropriate scaling. 
Thus S creates "squeezed cat states" out of coherent ones, taking a state localized at (q,p) to two 
that are localized at (2gmod l,p/2) and (2(/mod 1, (p + l)/2). Repeated action by S on an initial 
coherent state is illustrated in Fig. and exact revival occurs for the same reason that a deck of 
cards under the perfect riffle-shuffle reorders. 

Using the action of S we can construct a quantum baker's map. The action of choosing the left 
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FIG. 1: The correlation | (qp\S k \qn,po) | 2 as function of (q,p) for the case of N = 64, where \qp) is a toral 
coherent state localized at (q,p) On further applying S to the last figure produces the first as in this 
case S 6 is the identity. 



or right vertical partition is done by the projectors Pi and P 2 = In — Pi, where 




(5) 



The action of stretching and compression is implemented by S, which however produces an extra 
copy, shifted in momentum by one-half. Thus this is in the other horizontal partition that divides 
momentum into two equal halves. Thus we once again use projectors, now in momentum space 
to excise the extra copy and complete the action. The full quantum baker built around S is then 
written as: 



quantum baker's map since closer inspection shows that it is indeed very close to the usual baker's 

5 



B S = V2G~ N 



(PiG N SP 1 + P 2 G N SP 2 ) . 



(6) 



The factor of v2 is essential to restore unitarity after the projecting actions. This is not yet another 



map in Eq.lfT]). This is seen on rewriting Bg as: 

(- -) 

*N/2 

f\ 2 ' i> 



-((2 '4) n 



1 G N/2 

That the usual quantum baker's map is capable of generalizations, including arbitrary phases 
as boundary conditions and relative phases between the two blocks in the mixed representation 
is well-known |3j, though not all of these "decorated" baker's respect the symmetries of parity 
and time-reversal. The operator Bs however shows the explicit relationship between a quantum 
baker's map and the solvable operator S, whose action on the position basis is practically the 
doubling map restricted to the integers. It may be emphasized that even in Bs we are using anti- 
periodic boundary conditions, the phases of 1/4 and 3/4 in the Gjv/2 blocks (as well as the factor 
of i = v— l) is a direct consequence of the primitive structure in Eq. l(7j). That the operator obeys 
parity symmetry follows from the fact that Rn commutes with Gn and on verifying that 

^V/2^jv/2 = * B-N/2- (8) 

However it does not obey the time-reversal symmetry obeyed by the usual quantum baker's map. 
This follows from the preferential treatment of the position basis, in which S is a permutation, 
whereas in the momentum basis it is not. In the following we use S as an intermediate operator 
towards simplifying states of the usual baker's map B of Eq. lfl|). While doing so we will also 
compare the case of the operator Bs wherein there is a more explicit relationship; however a more 
detailed study of the spectra of Bs and related operators is itself postponed. 

The operator S is easily diagonalized. The case N = 2 K is particularly simple, as one sees from 
the cyclic shifting that S K = I/v, and therefore the possible eigenvalues are uj 1 where ujk = e 2m l K , 
and < I < K — 1. The complete set of eigenfunctions can be constructed based on the periodic 
orbits of the full binary left shift. When K is composite, an arbitrary isT-tuple may not produce (on 
action by S) an invariant subspace of full dimensionality K. Let the number of primitive periodic 
orbits of period n of the left shift map be denoted as p(n), this is the number of primitive binary 
n-tuples, where a primitive n-tuple is one that is not a repetition of a shorter string. If K has 
divisors di, da, ■ ■ ■ ,dM (including 1 and K), dimensionality of the invariant subspaces are di, and 
there are p(di) of them. In these subspaces the eigenfunctions maybe written as 

di-l 

rdi 



7^= Y. 4" V" a lL 2 ...a„ .. (9) 

m=0 

The corresponding eigenvalues u)~^. are p{di)-io\d degenerate. The number of primitive orbits 
is p{n) = J2k\n A* ( n /k) 2 k /n, where fi(n) is the Mobius function and the sum is over all the 
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FIG. 2: The quantum period function Tjv which is the multiplicative order of 2 modulo N — 1. 

divisors of n. A particularly simple case is when K is prime, as only d\ = 1 andcfo = K are the 
possible dimensions, and the states in the latter subspace have a degeneracy of (2 K — 2)/K = 
(N — 2)/K. Even when K is not prime the degeneracy increases in the same manner for large 
N. When N is not a power of 2, the matrix S has nontrivial spectral properties. SineeS \n) = 
|2*nmod(A r — 1)), there exists a time T/v such that S N = In- This must be the least integer 
such that2 T = 1 mod(A r — 1), the "quantum period function" Tv is then simply the multiplicative 
order defined above, Tv = ord,/v_i(2). This is not a simple function and its solution is equivalent 
to the difficult discrete logarithm problem, and thence to the task of factoring numbers.lt oscillates 
wildly with N, as seen in Fig. [21 going all the way from ln(iV)/ln(2) when TV is a power of 2 to 
(f>(N — 1) ~ (N — l)e _7 /ln(ln(iV — 1)), where 7 is the Euler constant. 

The eigenvalues are then T/y-th roots of unity, and one set of eigenfunctions are given by 



\<Pr) = ^ T jr\ W (^^) |2«mod(iV-l)), (10) 

where < r < T/v — 1. For certain iV the period T/v is maximal, that is T/v = ^(^V — 1) = N — 2. 
Naturally a necessary condition for this is that iV — 1 be prime. In this case apart from the 
eigenstates with unit eigenvalues, |0) and \N — 1), the others are exhaustively given by the above 
set. If T/v 7^ N — 2, other eigenfunctions can be found based on other subgroups. In general there 
is degeneracy and the states reside in some appropriate subspace. 

If we use the eigenstates of S as a basis for the eigenstates of the quantum baker's map, B, or 
Bs we find remarkable simplifications, as indeed these operators are "close" to each other. The 
crucial difference is that we can solve for the spectrum of S exactly. There are evident similarities 
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FIG. 3: The participation ratio in the position and S'-basis of the quantum baker's maps B (left) and Bs 
(right) when N = 198 (top) and N = 2222 (bottom). These cases are such that T N = N - 2. In all the 
figures the lower curve corresponds to the S-basis, while the upper one to the position basis. The states are 
arranged in the increasing order of the participation ratio in the S'-basis. 

of S to the well-studied quantum cat maps jl4|, w here there is a quantum period function that is 
wildly oscillatory, exactly solvable eigenstates Jl5|and so on. Here however the mathematics is far 
simpler, involving as it does a scalar multiplier (namely 2) rather than an integer 2x2 matrix. 

Let the eigenvectors of S be \(f> r ), we then evaluate the participation ratio (PR) l/Er K^|V')| )> 
which gives us (roughly) the number of S eigenstates needed to construct the vector \ip), here 
chosen to be one of the eigenstates of B. This is the PR in the S'-basis, while the PR in the 
position basis is similarly defined and indicates the derealization in position. For complex random 
states random matrix theory predicts a PR of N/2. In Fig. © we compare the participation ratio 
of the eigenstates of B and Bs in both the position and the S'-basis for a particular case, when 
the the S spectrum is largely non-degenerate. The PR in the position basis is halved to take into 
account the parity symmetry of the eigenstates, the S-basis already having this symmetry. We 
notice that the S-basis "simplifies" the states significantly as the PR is lesser by a factor of about 
five or more. 

We see from the figure that the S-basis simplifies states significantly more in the case of the 
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FIG. 4: The intensity of three cigcnstates of the quantum baker's map (N — 198) shown in both the position 
(left) and in the S'-basis (right). The states were chosen for their contrast in the two basis. 

operator B$, rather than the usual quantum baker's map. At the same time, large dips are seen 
for the eigenstates of B$ that are not visible for B, indicating perhaps that deviations from RMT 
(Random Matrix Theory) are larger in the case of the spectra of Bs- To illustrate the simplification, 
we show in Fig. Q three eigenstates of the usual quantum baker's map B, for the case N = 198 
that are considerably simplified in the S'-basis. 

We may improve upon the S'-basis by making it compliant with time-reversal symmetry. For 
instance, in the first state (say shown in Fig. (|4"|). the maximum overlap with an S-eigenstate 
\4> r ) is |(</v|V')| 2 = 0.34, while the (unnormalized) adapted state \4>' r ) = \(f> r ) + G~^\(j) r )* has an 
overlap of 0.37. This adapted state is such that Gj<r\<t>' r ) = \4>' T )* as required by time-reversal 
invariance of the quantum baker's map. An arbitrary phase between \<p r ) and Gj^~\<j) r )* was set 
as zero after numerically ascertaining that this was the optimal value. Note that the conjugation 
assume that the states are in the postion representation. 



We remark that this simplification fal 



basis for the case when A is a power of 2 |2|. In this case (for the operator B) the Thue-Morse states 
and many others simplified considerably more in the Hadamard basis, or after a Walsh-Hadamard 



s significantly short of that achieved by the Hadamard 
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transform; for instance in the case when N = 1024, after the transform the participation ratio 
of the Thue-Morse state was of the order of 2. While the Thue-Morse sequence (rather its finite 
truncations) is an eigenstate of S, the Hadamard transform itself commutes with S. Due to the 
degeneracy in the spectrum of S, it appears that the Hadamard transform represents a basis that 
is more optimal than that provided by the eigenvectors of S. The meaning of this commutation of 
S and H perhaps in terms of a classical symmetry is not clear to the author. 

Finally we remark on the statistical properties of the eigenstates, and on the "relative random- 
ness", in the sense of Kus and Zyczkowski 16], of S and the operators B and B$- The usual 
quantum baker's map eigenstates are nearly generic in the sense that they are close to those that 
are expected from RMT |17J, however there are also known and significant deviations, whose ori- 
gins may be number-theoretic (such as the multifractal scaling of eigenstates for the case when N 
are powers of 2 7j). We find, from results not presented here, that while the eigenstates in the 
position basis are much closer to the expected Porter-Thomas distribution, the eigenstates in the 
S'-basis are considerably deviated, as is to be expected. 

To quantitatively compare S and the baker's map operators B and B$, we study their relative 
randomness, or degree of noncommutativity, by means of the inner-product between the operator 
S and its image under B (or B$)- Thus define 

R 1 = \(S\BSB^)\/N, (11) 

where (X\Y) = Tt(X Y^). It is argued in 16] that this (and related quantities) are small, near zero, 
if the operators S and B are relatively random, whereas if they commute or anticommute R\ = 1. 
We show in Fig. (JSJ) this measure for both the operators B and B$ as a function of N. It is clear 
that the quantum baker's map B is significantly correlated to the operator S, as the inner-product 
R\ is around 0.4, and that the operator B$ is more so correlated, as the inner-product is around 
0.5. This is of course reflected in the fact that the eigenstates of B$ are more compressed in the 
/S'-basis. It is worthwhile remarking that powers of 2 do not appear to be special for the measure 
R\. Also the inner-products between S and the baker's map operators themselves behave similarly, 
as \{S\B)\/N « 0.63 while \{S\B S )\/N ra 0.70. 

In conclusion the exactly solvable operator S is a good "zeroth order" system for the quantum 
baker's map. This operator is somewhat similar to the semiquantum operators that are obtained 
on quantizing classical baker's after times larger than one [6J. However these operators usually 
have complicated spectra themselves. We can use S to build a quantum baker's map, which is very 
close to the usual baker's map, which in turn explains the close relationship between the solvable 



10 



0.54 
0.52 
0.5 
0.48 

of 

0.46 
0.44 
0.42 
0.4 

20 40 60 80 100 120 140 160 180 200 
N 

FIG. 5: The relative randomness measure R± as function of N, between the operators S and Bs (upper 
curve), and between S and B (lower curve). 

spectrum of S and that of the quantum baker's map. Using a relative randomness measure it has 
been shown that indeed the operator S is significantly correlated with the quantum baker's map. 
While pointing to the evident connection of S to the task of factoring numbers, it is tempting to 
speculate that the relationship between classically hard computations and their (probably faster) 
quantum algorithms has a deeper connection to the transition from classical to quantum chaos. 
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